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THE INTERMEDIATE-SCALE CLUSTERING OF LUMINOUS RED GALAXIES 
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ABSTRACT 

We report the intermediate-scale (0.3 to 40/i _1 Mpc) clustering of 35, 000 luminous early-type galaxies 
at redshifts 0.16 to 0.44 from the Sloan Digital Sky Survey. We present the redshift-space two-point 
correlation function £(s), the projected correlation function w p (r p ), and the deprojected real-space cor- 
relation function £(r), for approximately volume-limited samples. As expected, the galaxies are highly 
clustered, with the correlation length varying from 9.8 ± 0.2/i _1 Mpc to 11.2 ± 0.2ft," 1 Mpc, dependent 
on the specific luminosity range. For the —23.2 < M g < —21.2 sample, the inferred bias relative to that 
of galaxies is 1.84 ± 0.11 for 1ft -1 Mpc < r p < 10/i _1 Mpc, with yet stronger clustering on smaller 
scales. We detect luminosity-dependent bias within the sample but see no evidence for redshift evolution 
between z = 0.2 and z = 0.4. We find a clear indication for deviations from a power-law in the real-space 
correlation function, with a dip at ~ 2/i _1 Mpc scales and an upturn on smaller scales. The precision 
measurements of these clustering trends offer new avenues for the study of the formation and evolution 
of these massive galaxies. 

Subject headings: cosmology: observations — galaxies: clustering — galaxies: clusters: general — 

galaxies: distances and redshifts — galaxies: elliptical and lenticular, cD — galaxies: 
evolution — galaxies: statistics — large-scale structure of universe 



1. INTRODUCTION 

The clustering of galaxies provides a window not only 
to the formation of inhomogeneities in the early universe 
but also onto the physics of galaxy formation. Galaxies 
with different properties cluster differently (Hubble 1936; 
Zwicky et al. 1968; Davis & Geller 1976; Dressier 1980; 
Postman & Geller 1984; Hamilton 1988; White, Tully, & 
Davis 1988; Park et al. 1994; Loveday et al. 1995; Guzzo 
et al. 1997; Benoist et al. 1996; Willmer, da Costa & Pelle- 
grini 1998; Brown, Webster & Boyle 2000; Carlberg et al. 
2001; Norberg et al. 2001; Zehavi et al. 2002; Norberg et 
al. 2002; Budavari et al. 2003; Madgwick et al. 2003; Hogg 
et al. 2003; Zehavi et al. 2004b), and these trends can be 
connected to their small-scale environments, notably the 
masses of their host dark matter halos (e.g., Kaiser 1984; 
Bardeen et al. 1986; Mo & White 1996; Benson et al. 2000; 
Sheth, Mo, & Tormen 2001; Berlind et al. 2003). This path 
has been strengthened recently by the discovery of devi- 
ations from the canonical power-law correlation function 
on small scales (e.g., Zehavi et al. 2004a; Zheng 2004) and 
the ease of interpretation of these features by contempo- 
rary models of galaxy and halo clustering, in terms of the 
clustering of galaxies within single halos and the cluster- 
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mg between halos (Kauffmann, Nusser, & Steinmetz 1997; 
Jing, Mo, & Borner 1998; Kauffman et al. 1999; Benson 
et al. 2000; Ma & Fry 2000; Peacock & Smith 2000; Sel- 
jak 2000; Scoccimarro et al. 2001; Berlind & Weinberg 
2002; Berlind et al. 2003; Magliocchetti & Porciani 2003; 
Kravtsov et al. 2004; Zehavi et al. 2004a; Zheng 2004). 

The Sloan Digital Sky Survey (SDSS;York et al. 2000) 
was designed in scope and systematic control to permit 
the study of galaxy clustering over a wide range of scales 
and galaxy properties (e.g., Connolly et al. 2002; Zehavi et 
al. 2002; Budavari et al. 2003; Hogg et al. 2003; Tegmark 
et al. 2004; Zehavi et al. 2004b; all using the SDSS main 
galaxy sample). To improve the precision of clustering 
measurements on the largest scales, the SDSS provides 
a spectroscopic sample of luminous red galaxies (LRG). 
These galaxies reach a redshift of 0.5, thereby providing 
a sample of over lh~ 3 Gpc 3 (see Eisenstein et al. 2001). 
Thus far, over 50,000 spectra of LRGs have been acquired. 

In this paper, we will investigate the clustering of these 
luminous early-type galaxies on scales between 0.3 and 
40/i~ 1 Mpc. This stretches from the quasi- linear to the 
deeply non-linear regime. As massive early-type galax- 
ies are known to inhabit preferentially rich environments 
(e.g., Sandage 1972; Dressier 1980; Hoessel et al. 1980; 
Schneider et al. 1983; Postman & Geller 1984; Postman & 
Lauer 1995, M. Bernardi 2004, in preparation), this selec- 
tion should permit one to study the clustering and internal 
structure of massive halos. Models that differ in their asso- 
ciation of LRGs to cluster-sized halos or to the fraction in 
smaller halos will vary not only in their predicted correla- 
tion length, but also in the fine structure of the correlation 
functions. With the sample size available within the SDSS 
LRG sample, we expect to reach the precision necessary 
to perform such tests despite the rarity of massive galax- 
ies. We note that the rapid increase in the clustering of 
early-type galaxies at the highest luminosities (Hogg et al. 
2003; Zehavi et al. 2004b) implies that the connections be- 
tween the most massive galaxies and their environments is 
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notably different than even early-types. 

The outline of the paper is as follows. In § 2 we present 
the LRG sample. In § 3 we present the clustering mea- 
surements in redshift-space and in § 4 we show the inferred 
real-space clustering results. We conclude in § 5. Details 
of our sample modeling are given in the Appendix. 

Throughout the paper, all distances are comoving and 
quoted in /i _1 Mpc, where h = i/o/100 km s _1 Mpc -1 . 
For all distances and absolute magnitude we use a cosmol- 
ogy of O ro = 0.3 and A = 0.7 and adopt h = 1 to compute 
absolute magnitudes. 

2. DATA 

2.1. The SDSS LRG Sample 

The SDSS (York et al. 2000) is imaging 10 4 square de- 
grees away from the Galactic Plane in 5 passbands, it, 
g, r, i, and z (Fukugita et al. 1996; Gunn et al. 1998). 
Image processing (Lupton et al. 2001; Stoughton et al. 
2002; Pier et al. 2003) and calibration (Hogg et al. 2001; 
Smith et al. 2002) allow one to select galaxies, quasars, 
and stars for follow-up spectroscopy with twin fiber-fed 
double-spectrographs. The spectra cover 3800A to 9200A 
with a resolution of 1800. Targets are assigned to plug 
plates with a tiling algorithm that ensures nearly complete 
samples (Blanton et al. 2003a). An operational constraint 
of using fibers to obtain spectra is that no two fibers can 
be closer than 55" on the same plate. This constraint is 
partly alleviated by having roughly a third of the sky cov- 
ered by overlapping plates. 

Galaxy spectroscopic target selection proceeds by two 
algorithms. The primary sample (Strauss et al. 2002), 
referred to here as the MAIN sample, targets galaxies 
brighter than r < 17.77. The surface density of such 
galaxies is about 90 per square degree. The LRG algo- 
rithm (Eiscnstein et al. 2001) then selects ^12 additional 
galaxies per square degree, using color-magnitude cuts in 
g, r, and i to select galaxies to r < 19.5 that are likely to 
be luminous early- types at redshifts up to ~ 0.5. The selec- 
tion is extremely efficient, and the redshift success rate is 
very high. A few galaxies (3 per square degree at z > 0.15) 
matching the rest-frame color and luminosity properties of 
the LRGs arc extracted from the MAIN sample; we refer 
to this combined set as the LRG sample. In detail, there 
are two parts to the LRG algorithm, known as Cut I and 
Cut II and described in Eisenstein et al. (2001). 

We begin from a spectroscopic sample covering 3,836 
square degrees. The exact survey geometry is expressed in 
terms of spherical polygons and is known as lss_samplel4 
(M. Blanton 2004, in preparation). This set contains 55,000 
spectroscopic LRGs in the redshift range 0.15 < z < 0.55. 

2.2. Redshift and Magnitude Cuts 

The SDSS LRG sample is nearly volume-limited, but 
not precisely so. At z > 0.37, the flux limits of r < 19.2 
(Cut I) and r < 19.5 (Cut II) begin to move into the 
passively-evolving luminosity threshold. In this paper, we 
wish to analyze volume- limited samples, so as to study the 
clustering properties of well defined populations of galax- 
ies. We therefore define three subsamples in passively- 
evolved luminosity and restrict the redshift ranges to en- 
sure complete coverage. The subsamples are —23.2 < 
21.2 with 0.16 < z < 0.36, -23.2 < M g < -21.8 
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—21.2 sample. The shaded histogram is the distribution of the ac- 
tual data, and the solid continuous line is our model for the redshift 
distribution, described in Appendix A. The sample is close to a 
constant comoving volume for z < 0.36, although the fluctuations 
are reaching about 30% peak-to-peak (but one should note that the 
lowest redshifts, where the excess is, contain less volume than the 
redshift range would suggest). 
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Fig. 2.— As Figure 1, but for the -23.2 < M g < -21.8 sample. 
The sample is close to a constant comoving volume for z < 0.44. 



with 0.16 < z < 0.44, and -22.6 < M g < -21.6 with 
0.16 < z < 0.36. The first of these sets is picked to max- 
imize our use of the LRG spectroscopy for the innately 
volume-limited portion of the sample. The second is se- 
lected because the Cut II selection creates a knee in the 
number densities as a function of redshift that we can ex- 
ploit. The third is chosen to match the luminosity range 
of the -23 < M r < -22 volume-limited MAIN galaxy 
sample described in Zehavi et al. (2004b). We use here 
only the red subsample of the latter (as defined in Ze- 
havi et al. 2004b, with 0.10 < z < 0.23) to compare to 
the LRG clustering. These LRG and MAIN samples over- 
lap nearly completely in the redshift range in common, 
0.16 < z < 0.23. The basic information regarding the 
three LRG samples is summarized in Table 1, and their co- 
moving number density as a function of redshift is shown 
in Figures 1-3. 
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Table 1 



LRG Sample Statistics 





a 




Redshift 


Number 


Density b 


(M g ) c 




-23.2 < M g 


< 


-21.2 


0.16 < z < 0.36 


29298 


9.7 X 10" 5 


-21.63 


0.28 


-23.2 < M g 


< 


-21.8 


0.16 < z < 0.44 


12992 


2.4 x 10~ 5 


-22.01 


0.34 


-22.6 < M g 


< 


-21.6 


0.16 < z < 0.36 


14500 


4.8 x 10~ 5 


-21.84 


0.28 



a Rest-frame g-band absolute magnitudes, passively evolved to z = 0.3. 
b Average comoving densities are in units of h 3 Mpc — 3 . 
c Average rest-frame g-band absolute magnitude, M g 
d Average redshift 
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Fig. 3.— As Figure 1, but for the -22.6 < M g < -21.6 sample. 
The sample is close to a constant comoving volume for z < 0.36. 



The above luminosities have been fc-corrected and pas- 
sively evolved to rest-frame magnitudes at z = 0.3 (near 
the median redshift of the LRG sample) . We use the ob- 
served r-band to estimate rest- frame g, as this requires 
minimal ^-corrections at z = 0.3. We have used the "non- 
star-forming" model presented in Appendix B of Eisen- 
stein et al. (2001) but normalized to M g at z = 0.3. The 
model has relatively mild evolution, only about 1 magni- 
tude per unit redshift, compared to other measurements 
(Blanton et al. 2003b). Additional details of the samples' 
modeling are given in Appendix A. To the extent that our 
model is appropriate, our selections represent mass thresh- 
olds throughout the sample volume. However, as shown in 
Figure 1, the —23.2 < M g < —21.2 sample still has some 
redshift evolution in the number density. This is due to 
small fluctuations in the selection thresholds of the parent 
sample in luminosity and rest-frame color as a function 
of redshift. The other two samples, being safely more lu- 
minous than the LRG sample selection limits, are much 
closer to a constant comoving threshold (Figs. 2 and 3). 

3. REDSHIFT-SPACE CLUSTERING 

We calculate the LRG correlation function in redshift 
space as a function of the redshift-space separation s. To 



estimate the mean density and account for the complex 
survey geometry, we generate random catalogs, applying 
the radial and angular selection functions of the samples. 
The details of the radial and angular modeling are given 
in the Appendix. We typically use in each random catalog 
100-150 times the number of galaxies in the real sample, 
and we have verified that changing the random catalog 
makes negligible difference to the results. We estimate 
the correlation function using the Landy & Szalay (1993) 
estimator 

A DD -2DR + RR 

f = RR ' (1) 

where DD, DR and RR are the suitably normalized num- 
bers of weighted data-data, data-random and random-random 
pairs in each separation bin. We weight the galaxies (real 
and random) according to the angular and radial selec- 
tion functions. We use a simple weighting by the inverse 
of the selection function, as the samples we use are all 
approximately volume-limited, and we have verified that 
our results are insensitive to employing alternative weight- 
ing schemes. We also used the alternative £ estimators of 
Davis & Peebles (1983) and Hamilton (1993) and found 
no significant differences in the results. 

Here, and throughout the paper, we estimate statisti- 
cal errors on our measurements using jackknife resampling 
with 104 angular subsamples. Each subsample excludes 
roughly 37 square degrees (generally contiguous on the 
sky), which is about 90/i _1 Mpc comoving on a side at 
z = 0.3. The 2.5 degree SDSS stripes are 36/1" 1 Mpc co- 
moving at z = 0.3. Zehavi et al. (2004b) performed ex- 
tensive tests with mock catalogs to check the reliability of 
the jackknife error estimates over a similar range of sep- 
arations (see their Fig. 2). Their tests showed that the 
jackknife method is a robust way to estimate the error 
covariance matrix, especially for large volumes as probed 
here. 

Figure 4 shows the redshift-space correlation function, 
£(s), for the -23.2 < M g < -21.2 and -23.2 < M g < 
—21.8 LRG samples introduced in § 2.2, with errorbars 
obtained from the jackknife resampling. The small differ- 
ence in amplitude arises from the difference in the average 
luminosity of the galaxies in the samples (see Table 1), 
reflecting the known trend of stronger clustering with lu- 
minosity (e.g., Hamilton 1988; Park et al. 1994; Loveday 
et al. 1995; Benoist et al. 1996; Guzzo et al. 1997; Norberg 
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Fig. 5. — Redshift-space correlation function £(s) for the —22.6 < 
M g < —21.6 passively-evolved LRG sample and the comparable one 
from the MAIN galaxy sample. 



et al. 2001; Zehavi et al. 2002; Hogg et al. 2003; Zehavi et 
al. 2004b). The redshift-space correlation functions values 
are given in Table 2. 

Figure 5 shows £(s) for the -22.6 < M g < -21.6 LRG 
sample plotted together with £(s) obtained for the com- 
parable red —23 < M r < —22 MAIN galaxy subsample of 
Zehavi et al. (2004b), where we restrict the LRG sample to 
0.23 < z < 0.36, such that we are probing independent vol- 
umes. The LRG sample contains ~ 12400 galaxies, while 
the MAIN galaxy sample includes only ~2700. As is ob- 
vious from the plot, the agreement between the samples 
is excellent, and the LRG results thus extend in essence 
the MAIN galaxy clustering results to higher redshifts. 
The deviations at small separations are mainly due to shot 
noise effects arising from the small number of galaxies in 
the MAIN sample and are consistent within the errorbars. 
The numerical values of this LRG correlation function are 
also provided in Table 2. 



4. REAL-SPACE CLUSTERING 

4.1. Projected Correlation Function 

To separate effects of redshift distortions from spatial 
correlations, it is customary to estimate the correlation 
function on a two dimensional grid of pair separations 
parallel (tt) and perpendicular (r p ) to the line of sight, 
termed £(r p ,7r). One can then learn about the real-space 
correlation function by computing the projected correla- 
tion function 



■ p (r p ) = 2 dir£{r p ,ir). 



(2) 



In practice, we integrate up to TT max = 80ft. -1 Mpc, which 
is large enough to include most correlated pairs and gives 
a stable result. The omission of pairs at it > 80/i -1 Mpc 
likely causes an overestimate of w p (r p ) by l-2/i _1 Mpc, 
as the correlations on such scales are driven negative by 
redshift distortions. However, we have varied ir max from 
50-120/i _1 Mpc without significant change in w p . We also 
checked the robustness to binning in r p and in the integrated- 
over 7r direction, finding the results to be insensitive to 
either. 

The projected correlation function can in turn be related 
to the real-space correlation function, £(V), 



u p (r p ) = 2 rdr £(r)(r 2 - r/y 1 ' 2 



(3) 



(Davis & Peebles 1983). In particular, fitting a power-law 
to the w p (r p ) measurement allows us to infer the best-fit 
power law for £(r). 

Figure 6 shows the resulting projected correlation func- 
tion, w p (r p ), for the two inclusive LRG samples analyzed 
in this paper. Again, the differences in amplitude reflect 
the luminosity bias between the samples. Figure 7 shows 
the normalized (such that the diagonal is 1) jackknife er- 
ror covariance matrix of the w p (r p ) measurements for the 
—23.2 < M g < —21.2 sample. As is apparent, there is sig- 
nificant correlation between the measurements on different 
scales, but the auto-correlation along the diagonal is rela- 
tively strong with the cross-correlation falling off rapidly. 

Similar to the comparison to the MAIN galaxy sample 
results shown in § 3, Figure 8 compares the projected cor- 
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Table 2 



Correlation Functions Measurements 







23.2 < M g < - 


21.2 




23.2 < M g < - 


21.8 




22.6 < M g < - 


21.6 


separation 


CM 


w p (r p ) 


CM 


cm 


w p (r p ) 


CM 


CM 


w p (r p ) 


CM 


0.418 


87.4 (16.8) 


772.9 (52.7) 


675.7 (80.6) 


74.6 (44.0) 


1281 (243) 


1276 (399) 


39.6 (17.7) 


743.9 (86.8) 


484 (133) 


0.663 


37.7 (5.1) 


414.4 (26.0) 


210.3 (23.3) 


46.6 (17.5) 


575.2 (86.9) 


332.8 (84.8) 


30.3 (8.9) 


548.0 (66.8) 


316.9 (59.7) 


1.051 


25.1 (1.9) 


238.8 (13.4) 


66.9 (7.9) 


44.2 (8.3) 


286.8 (33.8) 


83.1 (21.8) 


30.9 (4.3) 


268.0 (29.0) 


75.5 (16.7) 


1.665 


18.5 (0.9) 


154.9 (7.5) 


24.9 (3.0) 


17.3 (2.3) 


178.5 (20.6) 


26.0 (7.8) 


19.4 (1.8) 


174.6 (14.1) 


30.8 (5.4) 


2.639 


11.0 (0.3) 


109.0 (5.1) 


11.7 (1.0) 


12.4 (1.0) 


132.4 (12.6) 


13.0 (3.2) 


12.1 (0.8) 


111.6 (9.0) 


10.1 (2.1) 


4.182 


6.32 (0.15) 


74.0 (3.8) 


5.04 (0.47) 


8.45 (0.45) 


95.5 (7.8) 


6.01 (1.11) 


7.35 (0.33) 


84.9 (7.1) 


5.58 (0.88) 


6.628 


2.99 (0.08) 


50.5 (2.6) 


2.32 (0.19) 


4.08 (0.19) 


70.8 (4.8) 


3.80 (0.50) 


3.52 (0.14) 


59.2 (4.1) 


2.62 (0.34) 


10.505 


1.28 (0.04) 


33.0 (2.2) 


1.04 (0.09) 


1.70 (0.08) 


38.8 (3.6) 


1.16 (0.19) 


1.52 (0.06) 


40.1 (3.1) 


1.25 (0.16) 


16.650 


0.54 (0.03) 


20.0 (1.8) 


0.43 (0.05) 


0.62 (0.03) 


25.2 (2.2) 


0.56 (0.09) 


0.59 (0.04) 


25.1 (2.5) 


0.57 (0.08) 


26.388 


0.19 (0.02) 


11.0 (1.4) 


0.14 (0.02) 


0.26 (0.02) 


13.2 (2.0) 


0.17 (0.06) 


0.23 (0.02) 


13.2 (1.8) 


0.20 (0.05) 



NOTES. — Measurements of the redshift-space correlation function, £(s), projected correlation function, w p (r p ), and real-space correlation 
function, §(r), for the three LRG samples discussed in the paper. Correlation functions are calculated for each sample over the range for which 
it is approximately volume-limited, denoted in Table 1. Comoving separations and w p (r p ) values are in /i -1 Mpc units. Redshift-space £(s) and 
real-space §(r) are dimensionless. The diagonal terms of the measurements error covariance matrices are given in parentheses. Our radial bins 
are logarithmically spaced with widths of 0.2 dex beginning at 10 — °' 49 . The separations listed in column 1 are the linear centers of the bins. 
Strictly speaking, the listed values of §(s) and w p (r p ) are the averages of these correlation functions over the annuli. However, for reasonable 
power-law interpolations, the values of w P (r p ) and £(r) are very nearly (<g; 1%) the values at the linear bin centers. The interpolated value of 
£(s) at the bin centers would be about 1.5% higher than the values in the table. 
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MAIN 0.10<z<0.23 
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Fig. 8. — Projected correlation function w p (r p ) for the —22.6 < 
M g < —21.6 LRG sample and the comparable one from the MAIN 
galaxy sample. 



Fig. 7. — Normalized error covariance matrix for the w p (r p ) mea- 
surement of the —23.2 < M g < —21.2 sample. The normalized 
covariance matrix is denned as Cij/(Cu ■ Cjj) 1 ^ 2 , where Cu are 
the elements of the error covariance matrix. Contour spacing is 0.1 
going from 1 on the diagonal (thick line) down to 0. The dashed 
lined denotes the 0.5 contour. Tickmarks denote the elements in 
the covariance matrix, and the labels denote the corresponding r p 
values. 



relation function of the analogous LRG and MAIN sam- 
ples. The small differences seen in the plot do not appear 
to be significant. A x 2 statistic of the difference, performed 
with the sum of the error covariance matrices of the two 
measurements, is 11.6 for the 10 degrees of freedom, con- 
sistent with cosmic variance. 

Now that we have demonstrated that the LRG sam- 



ple extends the MAIN galaxy sample, it is interesting to 
compare the amplitude of the LRG clustering to that of 
typical galaxies. Figure 9 shows such a comparison for 
our largest LRG sample (-23.2 < M g < -21.2) and the 
MAIN galaxy sample (a volume-limited sample with 
—21 < M r < —20 containing 5700 galaxies; Zehavi et al. 
2004b). Note that this is in contrast to the previous com- 
parisons (Fig. 5 and 8), where the LRG and red MAIN 
galaxy samples were chosen to match in luminosity and 
color. The quantities plotted are [w p (r p ) /w p ^ ld (r p )] 1 ^ 2 , 
where the fiducial u> p ^ ld (r p ) corresponds to a power-law 
correlation function £(r) = (r/5/i -1 Mpc) -1 ' 8 , and allow 
to infer the relative bias. For illustration purposes, we 
also plot this for a flat ACDM cosmology (with f2 m = 
0.3, h = 0.7, n = 1 and as — 0.9) projected correlation 
function computed from the nonlinear power spectrum of 
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Fig. 9. — [ w p( r p)/ w p^ ld ( r p)] 1 ^ 2 as a function of separation r p for 
the -23.2 < M g < -21.2 LRG sample (solid symbols and line) 
and for the L* MAIN galaxy sample (open symbols and short- 
dashed line; Zehavi et al. 2004b), allowing to infer their relative 
bias. Wpf td (r p ) is the projected correlation function corresponding 
to a fiducial power-law £(r) = (r/5/i _1 Mpc) - 1,8 . The long-dashed 
curve shows this relative quantity for a ACDM cosmology computed 
from the nonlinear power spectrum of Smith et al. (2003). 



Smith et al. (2003) (Z. Zheng, private communication). 
The matter correlation function is comparable in ampli- 
tude to the L.J, MAIN correlation function, but distinct 
in detail. For the LRG galaxies, this scaled quantity ap- 
pears roughly scale-invariant for 1 — 10/i _1 Mpc, with a 
notable upturn on smaller scales and a downturn on large 
scales. When fitting a constant bias factor between the 
two samples, taking into account the error covariance ma- 
trices, one obtains 6l_rg/^* = 1.84±0.11 when fitting over 
l/i -1 Mpc < r p < 10/i -1 Mpc. The scale dependence of 
the LRG inferred bias is in accord with the steeper correla- 
tion functions associated with red galaxies (e.g., Willmer, 
da Costa & Pellegrini 1998; Brown, Webster & Boyle 2000; 
Zehavi et al. 2002, 2004b). The downturn of the projected 
correlation function from a power-law on large separations 
is similar to that predicted by CDM models and to what 
is measured in the SDSS MAIN galaxy sample and in the 
2dF survey (Hawkins et al. 2003). 

Real-space correlation functions have been historically 
well described by power laws (Totsuji & Kihara 1969; Pee- 
bles 1974; Gott & Turner 1979; Davis & Peebles 1983; 
Fisher et al. 1994; Jing, Mo, & Borner 1998; Jing, Borner, 
& Suto 2002; Norberg et al. 2001; Zehavi et al. 2002), 
although recent precision measurements provide evidence 
for deviations from a power-law and a means of explaining 
them (e.g., Berlind et al. 2003; Magliocchetti & Porciani 
2003; Mailer et al. 2004; Zehavi et al. 2004a; Zheng 2004). 
Figure 10 shows power-law fits to our projected correla- 
tion functions. The inferred £(r) power-law fits are given 
in Table 3, while the w p (r p ) measurements themselves are 
provided in Table 2. Inspection of the values of the cor- 
relation length, ro, show clearly the trend with luminos- 
ity. The power-law slopes, 7, span the range 1.89 - 1.94. 
The x 2 /d.o.f values for the power-law fits are in the range 
2.3 — 3.9, indicating that a power-law is not a good fit. 
(The confidence level of a power-law fit is about 1.5% in 




: -23.2<M g <-21.8 0.16<z<0.44 

; -22.6<M g <-21.6 0.16<z<0.36 

- -33.3<M„<-21.2 0.16<z<0.36 
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Fig. 10. — Projected correlation function w p (r p ) for the three 
LRG samples discussed in the paper, plotted together with power- 
law fits fitted over the range 0.3 < r p < 30/i Mpc. 

Table 3 



£(r) Power Law Fits 



Mg 






ro 


7 


X 2 /d.o.f. 


-23.2 < Mg 


< 


-21.2 


9.80 ±0.20 


1.94 ±0.02 


3.9 


-23.2 < Mg 


< 


-21.8 


11.21 ±0.24 


1.92 ±0.03 


3.1 


-22.6 < Mg 


< 


-21.6 


10.59 ±0.29 


1.88 ±0.03 


2.3 



NOTES. — ro and 7 are obtained from a fit to w p (r p ) using the full 
error covariance matrix. The values of ro are given in h~ 1 Mpc units. 

the best case and less than 0.1% in the worst case.) 

The deviations from a power-law are clearly visible in 
Figure 11 where we divide the clustering measurements 
by a representative power-law w p (r p ) corresponding to a 
£(r) with ro — lOh^ 1 Mpc and 7 = 1.9. Similar deviations 
are also seen in the complementary analysis of the LRG 
samples by Eisenstein et al. (2004). These deviations ap- 
pear to be of a similar nature to the deviations detected in 
the MAIN galaxy samples (Zehavi et al. 2004a, b), which 
are naturally explained by contemporary models of galaxy 
clustering as the transition from a small-scale regime dom- 
inated by galaxy pairs in the same dark matter halo to a 
large-scale regime dominated by pairs of galaxies in sepa- 
rate halos. We delay to future work detailed halo modeling 
of this sort and interpretation of our measurements. There 
is a hint from the brightest sample of an increase at small 
scales (< l/i _1 Mpc) in the luminosity dependence of the 
bias, in agreement with the findings of Eisenstein et al. 
(2004). 



4.2. Luminosity and Redshift Dependences 

Figure 12 shows the redshift dependence of the —23.2 < 
M g < —21.8 results. One can see small deviations of 
the results corresponding to the different redshift ranges. 
For two independent redshift shells, we estimate the best- 
fitting multiplicative factor, a, between the two w p (r p ) 
measurements, taking into account the error covariance 
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Fig. 11. — Projected correlation function w p (r p ) for the LRG 
samples shown in Fig. 10, now all divided out by a fiducial power-law 
w p (r v ) corresponding to £(r) = (r/10/i — 1 Mpc) -1 ' 9 . The deviations 
from a power-law are clearly visible. 
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Fig. 12. — Redshift dependence of the w p (r p ) clustering results 
for the -23.2 < M g < -21.8 sample. 



matrices. This factor would be significantly different than 
one if redshift evolution was present and consistent with 
one otherwise. The multiplicative factor between the 0.16 < 
z < 0.23 and 0.23 < z < 0.36 results is a = 0.84±0.14 and 
between the 0.23 < z < 0.36 and 0.36 < z < 0.44 measure- 
ments itiset = 1.18±0.12. For our longest redshift base- 
line, we find a = 1.03 ± 0.17 between the 0.16 < z < 0.23 
and 0.36 < z < 0.44 measurements. Converting this to 
a limit on (1 + z) n , we find n = -0.2 ± 1.1. We thus 
conclude that while some variations between the different 
redshift shells are present, these are likely to reflect large- 
scale structure variations, and that no consistent trend of 
redshift evolution is detected. We note that for the w p (r p ) 
measurements for the comparable LRG and MAIN sam- 
ples shown in Figure 8, a — 1.01 ± 0.10, indicating clearly 
that no significant redshift evolution is present. 

We also checked the robustness of our results when cal- 
culating w p (r p ) separately for the two large disjoint areas 
of the northern Galactic sky covered in the current SDSS 
samples. The results from these two independent regions 



Fig. 13. — Luminosity dependence of w p (r p ) obtained for the 
three LRG samples for 0.16 < z < 0.36. 

are very similar and fully consistent, when calculated over 
the full redshift ranges of our samples. When looking at 
narrower redshift shells the differences tend to be a bit 
larger, reflecting the slight variations with redshift seen in 
Figure 12. For example, the tendency of w p (r p ) to have 
a slightly lower amplitude for the 0.16 < z < 0.23 shell is 
reproduced in the off-equatorial region, while w p (r p ) for 
the equatorial region is similar to that of the full volume. 
This supports our conclusion that these small deviations 
are sample variance effects reflecting the large-scale struc- 
ture fluctuations. 

Figure 13 shows the projected correlation function w p (r p ) 
obtained for the three LRG samples over an identical vol- 
umes (0.16 < z < 0.36). As mentioned previously, the 
small differences in clustering amplitude reflect the in- 
crease of clustering with luminosity. Again, we assess the 
significance of the increased clustering amplitude by esti- 
mating the best-fit multiplicative factor, a, between the 
measurements. Since these measurements are not fully in- 
dependent (they are obtained from the same volume and 
thus are susceptible to similar cosmic variance effects), we 
cannot treat their individual error covariance matrices as 
independent. Instead, we estimate the value of a from 
the mean and scatter of a obtained from the individual 
jackknife realizations. The resulting factor between the 
-23.2 < M g < 21.8 and -23.2 < M g < 21.2 measure- 
ments is a — 1.34 ± 0.08, more than a 4cr detection of 
luminosity bias. For the —22.6 < M g < —21.6 versus 
the -23.2 < M g < 21.2 measurements, a = 1.08 ± 0.05. 
It is clear that we detect a non-negligible luminosity bias 
among the different LRG samples. 

4.3. Real-Space Correlation Function 

It is possible to directly invert w p (r p ) to get £(r) in- 
dependent of the power-law assumption. This is done by 
recasting Equation 3 as 

f(r) = -i jH dr p d -^{r P 2 r^'\ (4) 

(e.g., Davis & Peebles 1983). We calculate the integral 
analytically by linearly interpolating between the binned 
w p (r p ) values, following Saunders et al. (1992). We note 
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-23.2<M <-81.i 

-22.6<M a <-21.6 

-23.2<M„<-21.2 



r (h 1 Mpc) 

Fig. 14. — Real-space correlation function £(r) for the three LRG 
samples. 



that this estimate is only accurate to a few percent level, 
due to the inaccuracy of the linear interpolation. Figure 14 
presents the real-space correlation function, obtained in 
this fashion, for the three LRG samples. The trends with 
luminosity and the hints of deviations from a power-law 
are noticeable here as well. The £(r) values for these sam- 
ples are given as well in Table 2. 

It is common to summarize the amplitude of the cor- 
relation function as the rms variation above Poisson in 
the counts of galaxies in R — 8/i _1 Mpc comoving radius 
spheres. The variance a 2 R can be calculated as 

4 = / d 3 n [ d 3 r 2 a\r2 - n\) (5) 

J\fi\<B, J\r 2 \<R 

where the integrals are over the interior of two spheres of 
radius R. This can be simplified to 

r2 / n„. o„,3^ 



'R 



/W>(3-| + f), 



(6) 



a useful formula that seems to have dropped out of the 
standard lore. Following Eisenstein (2003), we express this 
integral in terms of w p as 



n 2 1 

aR= W 



dr p r p w p (r p )g(r p /R) 



where g(x) is 
37r — 9x 



l 

1 

2tt 



■ X 
-28 



Vx 1 



+ a; 3 -9.T + 6sin _1 (2/a;) 



(7) 

for x < 2, 
for x > 2. 



The kernel g(x) is simpler than it looks: it starts positive, 
goes through zero at x w 1.28, and then returns to zero as 
— 1/irx 3 at large x. It is differentiable at x = 2. 

We compute as by using Equation 7 and assuming that 
r p w p (r p ) is constant in each bin. This yields as — 1.80 ± 
0.03 for the -23.2 < M g < -21.2 sample and a a = 2.06 ± 
0.05 for the -23.2 < M g < -21.8 sample, with the errors 
obtained from the jackknife subsamples. We stress that 
this is the real-space, non-linear erg, which should not be 
directly compared with the linear-regime as values that 
are typically quoted for CDM model normalization. 

To facilitate comparison with the cross-correlations be- 
tween LRGs and L* galaxies presented by Eisenstein et 



al. (2004), we also compute the following integral of the 
real-space correlation function: 

1 f°° 

' 4nr 2 dr £(r)W(r), (8) 



where 



V 



W{r) 



■ exp 



2a 2 



(9) 



where ao is a constant scale factor. The A statistic iso- 
lates about one octave of scale in the real-space correla- 
tion function. For power-law correlations of the observed 
slopes, A(ao) ~ £(1.75do). The A values are computed 
as a simple non-singular integral over w p (r p ). For the 
-23.2 < M g < -21.2 sample, we compute A = 66.5 ±3.5, 
15.6±0.7, 4.40±0.19, and 1.31±0.06 for a = 0.5k' 1 Mpc, 
lft,- 1 Mpc, 2/i~ 1 Mpc, and 4/i~ 1 Mpc proper distance at 
z = 0.3. These scales are chosen to match those in Eisen- 
stein et al. (2004). For the -23.2 < M g < -21.8 sample, 
the A values are 91.8 ± 9.6, 18.1 ± 1.8, 5.54 ± 0.43, and 
1.73 ± 0.12, respectively. 



5. CONCLUSIONS 

We have presented a statistical analysis of the intermediate- 
scale (0.3 to 40/i _1 Mpc) correlations of 35, 000 luminous 
red galaxies from the SDSS, using three nearly volume- 
limited subsamples. The size of the sample permits mea- 
surements of superb precision for these rare galaxies. We 
find clear deviations from power-law models in the pro- 
jected correlation functions. Relative to a power-law, there 
are excesses of clustering on sub-Mpc scales and at 5 — 
10 Mpc scales. These match qualitatively the deviations 
found in the SDSS MAIN sample galaxies and predicted 
by halo modeling (Zehavi et al. 2004a, b). 

The SDSS LRG sample reproduces the clustering re- 
sults obtained with the SDSS MAIN galaxy sample when 
matched appropriately in magnitude and color and thus 
provides an extension of the galaxy clustering analyses 
to higher redshifts, and a better signal-to- noise measure- 
ment of high luminosity galaxies. We find no evidence for 
redshift evolution of the correlation functions (for a fixed 
passively evolving magnitude range) out to z — 0.4. How- 
ever, the errors are such that dependences even as strong 
as (1 + z) ±2 cannot be excluded at more than 2a. 

All of the LRG samples are highly clustered, with corre- 
lation lengths around lOh Mpc comoving, roughly twice 
that of L* galaxies (e.g., Norberg et al. 2001; Zehavi et al. 
2002, 2004b). For the -23.2 < M g < -21.2 LRGs, the 
inferred bias relative to that of L* galaxies is 1.84 ±0.11 
on scales of 1 — lO/i" 1 Mpc. The bias is roughly scale in- 
variant on these scales and shows stronger clustering on 
smaller scales. 

We find that more luminous LRGs are yet more clus- 
tered; however, none of our samples reach the correlation 
levels of rich clusters (Bahcall & Soneira 1983; Nichol et 
al. 1992; Peacock & West 1992; Croft et al. 1997; Abadi 
et al. 1998; Lee & Park 1999; Collins et al. 2000; Gonza- 
lez, Zaritsky & Wechsler 2002; Bahcall et al. 2003). Note 
that all the latter estimate the cluster redshift-sp&ce cor- 
relation function, and thus should be compared to our So 
values inferred from Table 2 and not the ro values quoted 
in Table 3. The LRG clustering strengths and mean sep- 
arations d are comparable to those of the poorest clusters 
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mentioned in these works. Our measurements arc roughly 
consistent with the Bahcall et al. (2003) trend of increas- 
ing correlation length with mean separation, so = 2.6\/d 
(see also their Fig. 2). The LRG clustering strength we 
find is comparable as well to that of rich groups, which 
again have similar mean separations (Padilla et al. 2004; 
Yang et al. 2004). 

The SDSS LRG sample offers an enormous data set for 
the study of rare but important massive early-type galax- 
ies. The interplay of number density, clustering ampli- 
tude, correlation function shapes, redshift distortions, and 
higher-order correlations will provide a rich data set for 
the modeling of the relationship of these galaxies to their 
host halos and thereby to the evolution of the extreme end 
of the galaxy mass function. 
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APPENDIX A 
MODELING OF SELECTION FUNCTIONS 



A.l. Redshift Distribution 

To perform a 3-dimensional correlation analysis, one 
must have a model for the expected number of galaxies 
(in the absence of clustering) as a function of redshift at 
every point on the sky. Often this is done by computing a 
luminosity function and then integrating to a given depth. 
In the case of the LRG sample, this is hard to implement 
reliably because the luminosity function is so steep that 
minor variations in the k and e corrections are important 
and because the color selection does not select the identical 
region of rest-frame color-luminosity space at all redshifts. 

Instead, we construct an approximate model of the red- 
shift distribution based on models of the selection and then 
low-pass fit this model to the observed redshift distribu- 
tion. This removes power on very large radial separation 



scales (below k — 0.04/iMpc -1 ), but provides an excellent 
model for smaller scales. 

In detail, we predict the g — r and r — i colors of early- 
type galaxies as a function of redshift by convolving the 
observed average LRG spectrum (Eisenstein et al. 2003) 
with the SDSS system response. We then tweak that color- 
rcdshift relation onto the observed photometry with low- 
pass filtering. The spectral breaks of early-type galaxies 
create subtle features in the color-redshift relation that are 
currently difficult to resolve in the empirical photometry. 
At each redshift, we create a Monte Carlo set of colors 
by adding Gaussian noise to the mean color. We then 
find how luminous a galaxy would have to be to pass the 
LRG selection cut (this is a highly sensitive function of 
color) and integrate a hypothesized passively-evolved lu- 
minosity function, following Blanton et al. (2003b), to find 
the number of galaxies. Summing over the Monte Carlo 
points and multiplying by the cosmological volume gives 
the predicted real density of galaxies per redshift bin for 
a passively-evolving model population. This distribution 
is forced onto the observed redshift distribution using a 
boxcar smoothing length of 0.07 in redshift. Our resulting 
model is shown as the solid line in Figures 1,2, and 3. This 
more complex procedure is needed as, note, for example, 
the bump at z — 0.34 that is caused by the G band absorp- 
tion feature moving from the g to r filter. A simple boxcar 
smoothing of the redshift histogram does not preserve the 
height of this feature, but predicting the height from the 
average LRG spectrum produces an excellent match. 

A. 2. Angular Selection Function 

Selected galaxies can fail to receive spectra for three 
major reasons: 1) they might fall within the fiber collision 
radius 55" of another target and not get assigned a fiber, 
2) the plate to which they were assigned might not have 
been observed yet, and 3) they might have fallen outside 
of the boundary of any plate. With the adaptive tiling of 
the survey, it is very rare for an isolated object inside a 
plate boundary to fail to be placed on the plate. 

We weight the LRGs to account for fiber collisions by 
using a friends-of-friends grouping algorithm with a 55" 
linking length (provided within the lss_samplel4 pack- 
age) run on all of the galaxy and quasar targets. Within 
each collision group, we find the number of objects that 
did get assigned to a plate and divide by the total number. 
The inverse of this is the weighting assigned to the LRGs 
that did get assigned to plates. Note that to the extent 
that the selection of targets within the collision group was 
strictly random, this approach is the correct procedure for 
two-point clustering on separations larger than 55", even if 
the colliding objects are at very different redshifts. Imag- 
ine that there are two objects in close proximity and we 
wish to study the correlations with a third distant object 
(e.g., to count the pairs). In a perfect survey, one would 
have two pairs, each with a different redshift separation. In 
one version of the SDSS, one of the pairs would be counted, 
but with double weight; in another version, the other pair 
would be counted with double weight. The ensemble aver- 
age over many such examples is unbiased with respect to 
the result of the perfect survey. Only the pair correspond- 
ing to the two close objects is being mistreated (as it will 
never exist in this example) , but that pair only matters for 
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very small separation structure. In truth, the priorities of 
all targets are not quite equal, such that LRGs will always 
lose to quasar candidates, but the priority is equal to the 
dominant MAIN targets. The key angular separation of 
55" corresponds to 0.22/i _1 Mpc at z = 0.3, and we thus 
restrict our measurements to separations larger than that. 

The lss_samplel4 package provides the angular geome- 
try of the spectroscopic survey expressed in terms of spher- 
ical polygons. The geometry is complicated: the spectro- 
scopic plates are circular and overlap, while the imaging 
is in long strips on the sky and there are some overlap 
regions of certain plates that may not have been yet ob- 
served. The resulting spherical polygons track all these ef- 
fects and characterize the geometry in terms of "sectors" , 
each being a unique region of overlapping spectroscopic 
plates. In each sector, we count the number of possible 
targets (LRG, MAIN, and quasar), excluding those missed 
because of fiber collisions, and the number of these that 
did get a spectrum. We weight the galaxies by the ratio of 
these numbers. Typically, this ratio is very close to unity 
(mean of 1.007). It is large only in instances of regions 
covered by two plates where one of the plates has not yet 
been observed. To exclude these extreme cases, we cut 
our sample at a ratio of 1.67, which eliminates 550 LRG 
spectra. 

We also mask from the sample any regions that are close 
to bright (heavily saturated) stars. At present, we have not 
included any other masking. Notably, we are not explic- 
itly excluding regions because of high reddening or poor 
image quality; however, such regions are typically not tar- 
geted for spectroscopy and would already have been ex- 
cluded from the above geometry. We are not masking 
out small-scale imaging defects such as saturation trails 
or satellite tracks, nor avoiding the region around bright 
galaxies. These involve less than 1% of the survey area 
and have negligible effects on the strong clustering signal 
shown by LRGs. We treat the boundaries of the spectro- 
scopic plates as simple boundaries of the survey. In other 
words, we neglect any correlation between these real- world 
boundaries and the LRG distribution. Since the LRGs are 
a subdominant population with mild angular clustering, 
it is very unlikely that there is any important cross-talk 
between the LRG distribution and regions missed because 
of the adaptive tiling. 

A. 3. Color Calibration Issues 

A persistent worry with clustering of SDSS LRGs is the 
high sensitivity of the selection to errors in the photomet- 
ric calibration (Eiscnstcin et al. 2001). Errors of 1% in 
r — i create 8% fluctuations in the number density of the 
—23.2 < M g < —21.2 sample, although the other samples 
are more robust because their minimum luminosities are 
not close to the selection boundary and the selection vari- 
ations with absolute magnitude are three times less than 
those with r — i color. Obviously, errors in the photomet- 
ric zeropoints will be correlated on the sky, leaving long 
stripes of false density fluctuations on the sky. 

However, the SDSS has demonstrated 2% rms calibra- 
tions (Abazajian et al. 2004), which would produce only 
0.025 effects in the correlation function. Moreover, we find 
that this is an overestimate because the calibration er- 
rors are for the most part only correlated along the SDSS 



stripes, not between them. As the calibrations are applied 
by secondary patches that are only about 0.6 degrees wide 
(Stoughton et al. 2002), the calibration errors are sup- 
pressed on larger angular scales. We see no indication of 
any significant effect on the scales in this paper from cali- 
bration effects. 
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